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Abstract 



A search for narrow resonances at high mass in the dimuon and dielectron channels 
has been performed by the CMS experiment at the CERN LHC, using pp collision 
data recorded at^/s = 7 TeV. The event samples correspond to integrated luminosities 
of 40 pb^^ in the dimuon channel and 35 pb^^ in the dielectron channel. Heavy 
dilepton resonances are predicted in theoretical models with extra gauge bosons (Z') 
or as Kaluza-Klein graviton excitations (Gkk) in the Randall-Sundrum model. Upper 
limits on the inclusive cross section of Z'(Gkk) — ^ £^£^ relative to Z — > i^£^ are 
presented. These limits exclude at 95% confidence level a Z' with standard-model- 
like couplings below 1140GeV, the superstring-inspired Z'^^, below 887 GeV, and, for 
values of the coupling parameter fc/Mpi of 0.05 (0.1), Kaluza-Klein gravitons below 
855 (1079) GeV. 
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1 Introduction 

Many models of new physics predict the existence of narrow resonances, possibly at the TeV 
mass scale, that decay to a pair of charged leptons. This Letter describes a search for reso- 
nant signals that can be detected by the Compact Muon Solenoid (CMS) detector at the Large 
Hadron Collider (LHC) [IJ at CERN. The Sequential Standard Model Z^gj^ with standard- 
model-like couplings, the Z'^ predicted by grand unified theories [2], and Kaluza-Klein gravi- 
ton excitations arising in the Randall-Sundrum (RS) model of extra dimensions ||3l|4l were used 
as benchmarks. The RS model has two free parameters: the mass of the first graviton excitation 
and the coupling fc/ Mpi, where k is the curvature of the extra dimension and Mpi is the reduced 
effective Planck scale. Two values of the coupling parameter were considered: k/ Mpi = 0.05 
and 0.1. For a resonance mass of 1 TeV, the widths are 30, 6 and 3.5 (14) GeV for a Zgg^, Z'^, and 
Gkk with k/M-pi = 0.05 (0.1), respectively 

The results of searches for narrow Z' — > and Gkk ~^ £^i^ resonances in pp collisions 

at the Tevatron with over 5fb^^ of integrated luminosity at centre-of-mass energy of 1.96 TeV 
have previously been reported [Sl-ISl. Indirect constraints have been placed on the mass of the 
virtual Z' bosons by LEP-II experiments ll9]- IT2ll by examining the cross sections and angular 
distribution of dileptons and hadronic final states in e^e^ collisions. 

The results presented in this Letter were obtained from an analysis of data recorded in 2010, cor- 
responding to an integrated luminosity of 40 ± 4 pb^^ in the dimuon charmel, and 35 ± 4pb^^ 
in the dielectron channel, obtained from pp collisions at a centre-of-mass energy of 7 TeV. The 
total integrated luminosity used for the electron analysis is smaller than that for the muon 
analysis because of the tighter quality requirements imposed on the data. The search for reso- 
nances is based on a shape analysis of dilepton mass spectra, in order to be robust against un- 
certainties in the absolute background level. By examining the dilepton-mass spectrum from 
below the Z resonance to the highest mass events recorded, we obtain limits on the ratio of 
the production cross section times branching fraction for high-mass resonances to that of the Z. 
Using further input describing the dilepton mass dependence on effects of parton distribution 
functions (PDFs) and fc-factors, mass bounds are calculated for specific models. In addition, 
model-independent limit contours are determined in the two-parameter (c^^, Cj,) plane [131. Se- 
lected benchmark models for Z' production are illustrated in this plane, where where Cu and 
are model-dependent couplings of the Z' to up- and down-type quarks, respectively allowing 
lower bounds to be determined. 

2 The CMS Detector 

The central feature of the CMS l|T4l apparatus is a superconducting solenoid, of 6 m internal 
diameter, providing an axial field of 3.8 T. Within the field volume are the silicon pixel and strip 
trackers, the crystal electromagnetic calorimeter (ECAL) and the brass/ scintillator hadron cal- 
orimeter (HCAL). The endcap hadronic calorimeters are segmented in the z-direction. Muons 
are measured in gas-ionization detectors embedded in the steel return yoke. In addition to the 
barrel and endcap detectors, CMS has extensive forward calorimetry. 

CMS uses a right-handed coordinate system, with the origin at the nominal interaction point, 
the X-axis pointing to the centre of the LHC, the y-axis pointing up (perpendicular to the LHC 
plane), and the z-axis along the anticlockwise-beam direction. The polar angle, 0, is measured 
from the positive z-axis and the azimuthal angle, (p, is measured in the x-y plane. 

Muons are measured in the pseudorapidity range \ri\ < 2.4, with detection planes based on one 
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3 Electron and Muon Selection 



of three technologies: drift tubes in the barrel region, cathode strip chambers in the endcaps, 
and resistive plate chambers in the barrel and part of the endcaps. The inner tracker (silicon 
pixels and strips) detects charged particles within the pseudorapidity range \ < 2.5. 

The electromagnetic calorimeter consists of nearly 76 000 lead tungstate crystals which provide 
coverage in pseudorapidity \t]\ < 1.479 in the barrel region (EB, with crystal size At] = 0.0174 
and = 0.0174) and 1.479 < < 3.0 in the two endcap regions (EE, with somewhat larger 
crystals.). A preshower detector consisting of two planes of silicon sensors interleaved with a 
total of 3 Xq of lead is located in front of the EE. 

The first level (LI) of the CMS trigger system, composed of custom hardware processors, selects 
the most interesting events using information from the calorimeters and muon detectors. The 
High Level Trigger (HLT) processor farm further decreases the event rate employing the full 
event information, including the inner tracker. The muon selection algorithms in the HLT use 
information from the muon detectors and the silicon pixel and strip trackers. The electromag- 
netic (EM) selection algorithms use the energy deposits in the ECAL and HCAL; the electron 
selection in addition requires tracks matched to clusters. Events with muons or electromagnetic 
clusters with pj above LI and HLT thresholds are recorded. 

3 Electron and Muon Selection 

3.1 Triggers 

The events used in the dimuon channel analysis were collected using a single-muon trigger. 
The algorithm requires a muon candidate to be found in the muon detectors by the LI trigger. 
The candidate track is then matched to a silicon tracker track, forming an HLT muon. The HLT 
muon is required to have px > 9 to 15 GeV, depending on the running period. 

A double EM cluster trigger was used to select the events for the dielectron channel. ECAL 
clusters are formed by summing energy deposits in crystals surrounding a "seed" that is lo- 
cally the highest-energy crystal. The clustering algorithm takes into account the emission of 
bremsstrahlung. This trigger requires two clusters with the ECAL transverse energy Ej above 
a threshold of 17 to 22 GeV, depending on the running period. For each of these clusters, the 
ratio H/E, where E is the energy of the ECAL cluster and H is the energy in the HCAL cells 
situated behind it, is required to be less than 15%. At least one of these clusters must have been 
associated with an energy deposit identified by the LI trigger. 

3.2 Lepton Reconstruction 

The reconstruction, identification, and calibration of muons and electrons follow standard CMS 
methods [15J. Combinations of test beam, cosmic ray muons, and data from proton collisions 
have been used to calibrate the relevant detector systems for both muons and electrons. 

Muons are reconstructed independently as tracks in both the muon detectors and the silicon 
tracker |T6| . The two tracks can be matched and fitted simultaneously to form a "global muon". 
Both muons in the event must be identified as global muons, with at least 10 hits in the silicon 
tracker and with pj > 20 GeV. All muon candidates that satisfy these criteria are classified 
as "loose" muons. At least one of the two muons in each event must be further classified as a 
"tight" muon by passing the following additional requirements: a transverse impact parameter 
with respect to the collision point less than 0.2 cm; a X'^ per degree of freedom less than 10 for 
the global track fit; at least one hit in the pixel detector; hits from the muon tracking system in 
at least two muon stations on the track; and correspondence with the single-muon trigger. 



3.2 Lepton Reconstruction 
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Electrons are reconstructed by associating a cluster in the ECAL with a track in the tracker l|T7l . 
Track reconstruction, which is specific to electrons to account for bremsstrahlung emission, is 
seeded from the clusters in the ECAL, first using the cluster position and energy to search for 
compatible hits in the pixel detector, and then using these hits as seeds to reconstruct a track 
in the silicon tracker. A minimum of five hits is required on each track. Electron candidates 
are required to be within the barrel or endcap acceptance regions, with pseudorapidities of 
< 1.442 and 1.560 < < 2.5, respectively. A candidate electron is required to deposit 
most of its energy in the ECAL and relatively little in the HCAL (H/E < 5%). The transverse 
shape of the energy deposit is required to be consistent with that expected for an electron, 
and the associated track must be well-matched in // and (p. Electron candidates must have 
£t > 25 GeV. 

In order to suppress misidentified leptons from jets and non-prompt muons from hadron de- 
cays, both lepton selections impose isolation requirements. Candidate leptons are required to 
be isolated within a narrow cone of radius AR = a/ (A;/)^ + (A^)^ = 0.3, centred on the lep- 
ton. Muon isolation requires that the sum of the px of all tracks within the cone, excluding the 
muon, is less than 10% of the pj of the muon. For electrons, the sum of the pj of the tracks, 
excluding the tracks within an inner cone of AR = 0.04, is required to be less than 7 GeV for can- 
didates reconstructed within the barrel acceptance and 15 GeV within the endcap acceptance. 
The calorimeter isolation requirement for electron candidates within the barrel acceptance is 
that, excluding the Ej of the candidate, the sum of the Ej resulting from deposits in the ECAL 
and the HCAL within a cone of AR = 0.3 be less than 0.03£t + 2 GeV. For candidates within 
the endcap acceptance, the segmentation of the HCAL in the z-direction is exploited. For can- 
didates with Ej below 50 GeV (above 50 GeV), the isolation energy is required to be less than 
2.5 GeV (0.03 (£t ~ 50) + 2.5 GeV), where Ej is determined using the ECAL and the first layer 
of the segmented HCAL. The Ej in the second layer of the HCAL is required to be less than 
0.5 GeV. These requirements ensure that the candidate electrons are well-measured and have 
minimal contamination from jets. 

The performance of the detector systems for the data sample presented in this paper is es- 
tablished using measurements of standard model (SM) W and Z processes with leptonic final 
states 1 15 1 and using traversing cosmic ray muons [18|. 

Muon momentum resolution varies from 1% at momenta of a few tens of GeV to 10% at mo- 
menta of several hundred GeV, as verified with measurements made with cosmic rays. The 
alignment of the muon and inner tracking systems is important for obtaining the best momen- 
tum resolution, and hence mass resolution, particularly at the high masses relevant to the Z' 
search. An additional contribution to the momentum resolution arises from the presence of 
distortion modes in the tracker geometry that are not completely constrained by the alignment 
procedures. The dimuon mass resolution is estimated to have an rms of 5.8% at 500 GeV and 
9.6% at 1 TeV. 

The ECAL has an ultimate energy resolution of better than 0.5% for unconverted photons with 
transverse energies above 100 GeV. The ECAL energy resolution obtained thus far is on average 
1.0% for the barrel and 4.0% for the endcaps. The mass resolution is estimated to be 1.3% at 
500 GeV and 1.1% at 1 TeV. Electrons from W and Z bosons were used to calibrate ECAL energy 
measurements. For both muons and electrons, the energy scale is set using the Z mass peak, 
except for electrons in the barrel section of the ECAL, where the energy scale is set using neutral 
pions, and then checked using the Z mass peak. The ECAL energy scale uncertainty is 1% in 
the barrel and 3% in the endcaps. 
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4 Event Samples and Selection 



3.3 Efficiency Estimation 

The efficiency for identifying and reconstructing lepton candidates is measured with the tag- 
and-probe method |T5| . A tag lepton is estabhshed by applying tight cuts to one lepton candi- 
date; the other candidate is used as a probe. A large sample of high-purity probes is obtained 
by requiring that the tag-and-probe pair have an invariant mass consistent with the Z boson 
mass (80 < niff < 100 GeV). The factors contributing to the overall efficiency are measured in 
the data. They are: the trigger efficiency, the reconstruction efficiency in the silicon tracker, the 
electron clustering efficiency, and the lepton reconstruction and identification efficiency. All 
efficiencies and scale factors quoted below are computed using events in the Z mass region. 

The trigger efficiencies are defined relative to the full offline lepton requirements. For the 
dimuon events, the efficiency of the single muon trigger with respect to loose muons is mea- 
sured to be 89% ± 2% IITSl . The overall efficiency, defined with respect to particles within the 
physical acceptance of the detector, for loose (tight) muons is measured to be 94.1% ± 1.0% 
(81.2% ± 1.0%). Within the statistical precision allowed by the current data sample, the dimuon 
efficiency is constant as a function of pj above 20 GeV, as is the ratio of the efficiency in the data 
to that in the Monte Carlo (MC) of 0.977 ± 0.004. For dielectron events, the double EM cluster 
trigger is 100% efficient (99% during the early running period). The total electron identification 
efficiency is 90.1% ± 0.5% (barrel) and 87.2% ± 0.9% (endcap). The ratio of the electron effi- 
ciency measured from the data to that determined from MC simulation at the Z resonance is 
0.979 ± 0.006 (EB) and 0.993 ± 0.011 (EE). To determine the efficiency applicable to high-energy 
electrons in the data sample, this correction factor is applied to the efficiency found using MC 
simulation. The efficiency of electron identification increases as a function of the electron trans- 
verse energy until it becomes flat beyond an Ej value of about 45 GeV. Between 30 and 45 GeV 
it increases by about 5%. 



4 Event Samples and Selection 

Simulated event samples for the signal and associated backgrounds were generated with the 
PYTHIA V6.422 mi MC event generator, and with MadGraph Jl^ and powheg vl.l 
123] interfaced with the pythia parton-shower generator using the CTEQ6L1 [24] PDF set. 
The response of the detector was simulated in detail using GEANT4 |.25il . These samples were 
further processed through the trigger emulation and event reconstruction chain of the CMS 
experiment. 

For both dimuon and dielectron final states, two isolated same flavour leptons that pass the 



lepton identification criteria described in Section 3.2 are required. The two charges are required 
to have opposite sign in the case of dimuons (for which a charge misassignment implies a large 
momentum measurement error), but not in the case of dielectrons (for which charge assignment 
is decoupled from the ECAL-based energy measurement). An opposite-charge requirement for 
dielectrons would lead to a loss of signal efficiency of a few percent. 

Of the two muons selected, one is required to satisfy the "tight" criteria. The electron sample 
requires at least one electron candidate in the barrel because events with both electrons in the 
endcaps will have a lower signal-to-background ratio. For both channels, each event is required 
to have a reconstructed vertex with at least four associated tracks, located less than 2 cm from 
the centre of the detector in the direction transverse to the beam and less than 24 cm in the di- 
rection along the beam. This requirement provides protection against cosmic rays. Additional 
suppression of cosmic ray muons is obtained by requiring the three-dimensional opening angle 
between the two muons to be smaller than n — 0.02 radians. 
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5 Backgrounds 

The most prominent SM process that contributes to the dimuon and dielectron invariant mass 
spectra is Drell-Yan production (Z/7*); there are also contributions from the tt, tW, WW, and 
Z — ?► TT channels. In addition, jets may be misidentified as leptons and contribute to the 
dilepton invariant mass spectrum through multi-jet and vector boson + jet final states. 

5.1 Z/7* Backgrounds 

The shape of the dilepton invariant mass spectrum is obtained from Drell-Yan production us- 
ing a MC simulation based on the pythia event generator. The simulated spectrum at the 
invariant mass peak of the Z boson is normalized to the data. The dimuon analysis uses the 
data events in the Z mass interval of 60-120 GeV; the dielectron analysis uses data events in the 
narrower interval of 80-100 GeV in order to obtain a comparably small background contamina- 
tion. 

A contribution to the uncertainty attributed to the extrapolation of the event yield and the 
shape of the Drell-Yan background to high invariant masses arises from higher order QCD 
corrections. The next-to-next-to-leading order (NNLO) /c-factor is computed using FEWZz 
vl.X |26|, with PYTHIA V6.409 and CTEQ6.1 PDF |j2Z| as a baseline. It is found that the vari- 
ation of the A;-factor with mass does not exceed 4% where the main difference arises from the 
comparison of PYTHIA and FEWZz calculations. A further source of uncertainty arises from 
the PDFs. The lhaglue 128J interface to the lhapdf-5.3.1 [29| library is used to evaluate these 
uncertainties, using the error PDFs from the CTEQ6.1 and the MRST2006nnlo ||30| uncertainty 
eigenvector sets. The uncertainty on the ratio of the background in the high-mass region to that 
in the region of the Z peak is below 4% for both PDF sets and masses below 1 TeV. Combining 
the higher order QCD and PDF uncertainties in quadrature, the resulting uncertainty in the 
number of events normalized to those expected at the Z peak is about 5.7% for masses between 
200 GeV and 1 TeV. 

5.2 Other Backgrounds with Prompt Lepton Pairs 

The dominant non-Drell-Yan electroweak contribution to high niu masses is tt; in addition 
there are contributions from tW and diboson production. In the Z peak region, Z ^ tt decays 
also contribute. All these processes are flavour symmetric and produce twice as many e}i pairs 
as ee or }i}i pairs. The invariant mass spectrum from e^f/^ events is expected to have the same 
shape as that of same flavour events but without significant contamination from Drell- 
Yan production. 

Figure [l] shows the observed e^jv^ dilepton invariant mass spectrum from a dataset corre- 
sponding to 35 pb^^, overlaid on the prediction from simulated background processes. This 
spectrum was obtained using the same single-muon trigger as in the dimuon analysis and by 
requiring oppositely charged leptons of different flavour. Using an electron trigger, a very sim- 
ilar spectrum is produced. Differences in the geometric acceptances and efficiencies result in 
the predicted ratios of /^^/^^ and ee to e^f/^ being approximately 0.64 and 0.50, respectively. In 
the data, shown in Fig.[l| there are 32 (7) e^ji^ events with invariant mass above 120 (200) GeV. 
This yields an expectation of about 20 (4) dimuon events and 16 (4) dielectron events. A direct 
estimate from MC simulations of the processes involved predicts 20.1 ± 3.6 (5.3 ± 1.0) dimuon 
events and 13.2 ± 2.4 (3.5 ± 0.6) dielectron events. The uncertainty includes both statistical 
and systematic contributions, and is dominated by the theoretical uncertainty of 15% on the tt 
production cross section fSTllS^I. The agreement between the observed and predicted distribu- 
tions provides a validation of the estimated contributions from the backgrounds from prompt 
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5 Backgrounds 




Figure 1: The observed opposite-sign e^fi"^ dilepton invariant mass spectrum (data points). 
The uncertainties on the data points (statistical only) represent 68% confidence intervals for the 
Poisson means. Filled histograms show contributions to the spectrum from tt, other sources of 
prompt leptons (tW, diboson production, Z tt), and the multi-jet background (from Monte 
Carlo simulation). 

leptons obtained using MC simulations. 

5.3 Events with Misidentified and Non-Prompt Leptons 

A further source of background arises when objects are falsely identified as prompt leptons. 
The misidentification of jets as leptons, the principle source of such backgrounds, is more likely 
to occur for electrons than for muons. 

Backgrounds arising from jets that are misidentified as electrons include W — )■ ev + jet events 
with one jet misidentified as a prompt electron, and also multi-jet events with two jets misiden- 
tified as prompt electrons. A prescaled single EM cluster trigger is used for collecting a sample 
of events to determine the rate of jets misreconstructed as electrons and to estimate the back- 
grounds from misidentified electrons. The events in this sample are required to have no more 
than one reconstructed electron, and missing transverse energy of less than 20 GeV, to sup- 
press the contribution from Z and W events respectively. The probability for an EM cluster 
with H/E < 5% to be reconstructed as an electron is determined in bins of Ej and r] from a 
data sample dominated by multi-jet events and is used to weight appropriately events which 
have two such clusters passing the double EM trigger. The estimated background contribution 
to the dielectron mass spectrum due to misidentified jets is 8.6±3.4 (2.1 ±0.8) for mee > 120 
(200) GeV. 

In order to estimate the residual contribution from background events with at least one non- 
prompt or misidentified muon, events are selected from the data sample with single muons that 
pass all selection cuts except the isolation requirement. A map is created, showing the isolation 
probability for these muons as a function of px and rj. This probability map is corrected for 
the expected contribution from events with single prompt muons from tt and W decays and 
for the observed correlation between the probabilities for two muons in the same event. The 
probability map is used to predict the number of background events with two isolated muons 



5.4 Cosmic Ray Muon Backgrounds 
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based on the sample of events that have two non-isolated muons. This procedure, which has 
been validated using simulated events, predicts a mean background to m^^, > 120 (200) GeV of 
0.8 ±0.2 (0.2 ±0.1). 

As the signal sample includes the requirement that the muons in the pair have opposite electric 
charge, a further cross-check of the estimate is performed using events with two isolated muons 
of the same charge. There are no events with same-charge muon pairs and m^,^, > 120 GeV, a 
result which is statistically compatible with both the figure of 1.6 ± 0.3 events predicted from 
SM process using MC simulation and the figure of 0.4 ± 0.1 events obtained using methods 
based on data. 

5.4 Cosmic Ray Muon Backgrounds 

The ji^ji^ data sample is susceptible to contamination from traversing cosmic ray muons, 
which may be misreconstructed as a pair of oppositely charged, high-momentum muons. 
Cosmic ray events can be removed from the data sample because of their distinct topology 
(collinearity of two tracks associated with the same muon), and their uniform distribution 
of impact parameters with respect to the collision vertex. The residual mean expected back- 
ground from cosmic ray muons is measured using sidebands to be less than 0.1 events with 
m^^ > 120 GeV. 

6 Dilepton Invariant Mass Spectra 

The measured dimuon and dielectron invariant mass spectra are displayed in Figs. |2jleft) and 
(right) respectively, along with the expected signal from Zggj^ with a mass of 750 GeV. In the 
dimuon sample, the highest invariant mass event has m^,^ = 463 GeV, with the pj of the two 
muons measured to be 258 and 185 GeV. The highest invariant mass event in the dielectron 
sample has nigg = 419 GeV, with the electron candidates having Ej of 125 and 84 GeV. 

The expectations from the various background sources, Z/7*, tt, other sources of prompt lep- 
tons (tW, diboson production, Z — > tt) and multi-jet events are also overlaid in Fig.|2j For the 
dielectron sample, the multi-jet background estimate was obtained directly from the data. The 
prediction for Drell-Yan production of Z/ 7* is normalized to the observed Z — > signal. All 
other MC predictions are normalized to the expected cross sections. Figures |3]^left) and (right) 
show the corresponding cumulative distributions of the spectra for the dimuon and dielectron 
samples. Good agreement is observed between data and the expectation from SM processes 
over the mass region above the Z peak. 

Searches for narrow resonances at the Tevatron B |8| have placed lower limits in the mass 
range 600 GeV to 1000 GeV. The region with dilepton masses 120 GeV < niff < 200 GeV is 
part of the region for which resonances have been excluded by previous experiments, and thus 
should be dominated by SM processes. The observed good agreement between the data and 
the prediction in this control region gives confidence that the SM expectations and the detector 
performance are well understood. 

In the Z peak mass region defined as 60 < nifj- < 120 GeV, the number of dimuon and dielec- 
tron candidates are 16 515 and 8 768 respectively, with very small backgrounds. The difference 
in the electron and muon numbers is due to the higher Ej cut in the electron analysis and lower 
electron identification efficiencies at these energies. The expected yields in the control region 
(120-200 GeV) and high invariant mass regions (> 200 GeV) are listed in Table [l] The agree- 
ment between the observed data and expectations, while not used in the shape-based analysis, 
is good. 
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6 Dilepton Invariant Mass Spectra 



CMS \/s = 7TeV Ldt = 40pb' CMS \/s = 7TeV L dt = 35 pb' 




m{\xY) [GeV] m(ee) [GeV] 



Figure 2: Invariant mass spectrum of (left) and ee (right) events. The points with error 

bars represent the data. The uncertainties on the data points (statistical only) represent 68% 
confidence intervals for the Poisson means. The filled histograms represent the expectations 
from SM processes: Z/7*, tt, other sources of prompt leptons (tW, diboson production, Z 
tt), and the multi-jet backgrounds. The open histogram shows the signal expected for a Zggjy[ 
with a mass of 750 GeV. 




Figure 3: Cumulative distribution of invariant mass spectrum of (left) and ee (right) 

events. The points with error bars represent the data, and the filled histogram represents the 
expectations from SM processes. 
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Table 1: Number of dilepton events with invariant mass in the control region 
120 < < 200 GeV and the search region m^^ > 200 GeV. The expected number of 

Z' events is given within ranges of 328 GeV and 120 GeV for the dimuon sample and the di- 
electron sample respectively, centred on 750 GeV. The total background is the sum of the SM 
processes listed. The MC yields are normalized to the expected cross sections. Uncertainties 
include both statistical and systematic components added in quadrature. 



Source 




Number of events 






Dimuon sample 


Dielectron sample 




(120 -200) GeV 


>200GeV 


(120 - 200) GeV 


>200GeV 


CMS data 


227 


35 


109 


26 


ZssM(750GeV) 




15.0 ±1.9 




8.7 ± 1.1 


Total background 


204 ± 23 


36.3 ±4.3 


120 ± 14 


24.4 ± 3.0 


Z/7* 


187 ± 23 


30.2 ± 3.6 


104 ± 14 


18.8 ±2.3 


tt 


12.3 ± 2.3 


4.2 ±0.8 


7.6 ± 1.4 


2.7 ±0.5 


Other prompt leptons 


4.4 ± 0.5 


1.7 ±0.2 


2.1 ±0.2 


0.8 ±0.1 


Multi-jet events 


0.6 ±0.2 


0.2 ±0.1 


6.5 ±2.6 


2.1 ±0.8 



7 Limits on the Production Cross Section 

The observed invariant mass spectrum agrees with expectations based on standard model pro- 
cesses, therefore limits are set on the possible contributions from a narrow heavy resonance. 
The parameter of interest is the ratio of the products of cross sections and branching fractions: 

_ c7(pp^Z^±X^^f±X) 
" a{pp ^ Z ± X ^ ± X) ■ ^ ' 

By focusing on the ratio Ra, we eliminate the uncertainty in the integrated luminosity, reduce 
the dependence on experimental acceptance, trigger, and offline efficiencies, and generally ob- 
tain a more robust result. 

For statistical inference about Ra, we first estimate the Poisson mean jiz of the number of 
Z — )■ events in the sample by counting the number of events in the Z peak mass region 
and correcting for a small (~ 0.4%) background contamination (determined with MC simula- 
tion). The uncertainty on jiz is about 1% (almost all statistical) and contributes negligibly to the 
uncertainty on R^. 

We then construct an extended unbinned likelihood function for the spectrum of ££ invariant 
mass values m above 200 GeV, based on a sum of analytic probability density functions (pdfs) 
for the signal and background shapes. 

The pdf /s {tn\T, M, w) for the resonance signal is a Breit-Wigner of width T and mass M convo- 
luted with a Gaussian resolution function of width w (section [3.2| |. The width T is taken to be 
that of the Zggj^ (about 3%); as noted below, the high-mass limits are insensitive to this width. 
The Poisson mean of the yield is /^s = Ra ■ }iz ■ Re, where Re is the ratio of selection efficiency 
times detector acceptance for Z' decay to that of Z decay; }Ib denotes the Poisson mean of the 
total background yield. A background pdf /b was chosen and its shape parameters fixed by fit- 
ting to the simulated Drell-Yan spectrum in the mass range 200 < mg£ < 2000 GeV. Two func- 
tional forms for the dependence of /b on shape parameters a and k were tried: fB{m\<x,K) ~ 
exp(— am*^) and ~ exp(— am)m^'^ Both yielded good fits and consistent results for both the 
dimuon and dielectron spectra. For definiteness, this Letter presents results obtained with the 
latter form. 
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7 Limits on the Production Cross Section 



The extended likelihood L is then 

where m denotes the dataset in which the observables are the invariant mass values of the 
lepton pairs, nii; N denotes the total number of events observed above 200 GeV; and }i = }is + 
jis is the mean of the Poisson distribution from which N is an observation. 

Starting from Eqn. |2| confidence /credible intervals are computed using more than one ap- 
proach, both frequentist (using profile likelihood ratios) and Bayesian (multiplying C by prior 
pdfs including a uniform prior for the signal mean). With no candidate events in the region of 
small expected background above 465 GeV, the result is insensitive to the statistical technique, 
and also with respect to the width of the Z' and to changes in systematic uncertainties and their 
functional forms, taken to be log-normal distributions with fractional uncertainties. 

For Re, we assign an uncertainty of 8% for the dielectron channel and 3% for the dimuon 
channel. These values reflect our current understanding of the detector acceptance and recon- 
struction efficiency turn-on at low mass (including PDF uncertainties and mass-dependence 
of /c-factors), as well as the corresponding values at high mass, where cosmic ray muons are 
available to study muon performance but not electron performance. The uncertainty in the 
mass scale affects only the mass region below 500 GeV where there are events in both channels 
extrapolating from the well-calibrated observed resonances. For the dielectron channel, it is 
set to 1% based on linearity studies. For the dimuon channel, it is set to zero, as a sensitivity 
study showed negligible change in the results up to the maximum misalignment consistent 
with alignment studies (corresponding to several percent change in momentum scale). The 
acceptance for Gkk (spin 2) is higher than for Z' (spin 1) by less than 8% over the mass range 
0.75-1.1 TeV This was conservatively neglected when calculating the limits. 

In the frequentist calculation, the mean background level /^b is the maximum likelihood esti- 
mate; in the fully Bayesian calculation a prior must be assigned to the mean background level, 
but the result is insensitive to reasonable choices (i.e., for which the likelihood dominates the 
prior). 

The upper limits on Ra (Eqn. [T]) from the various approaches are similar, and we report the 
Bayesian result (implemented with Markov Chain Monte Carlo in RooStats p3| ) for definite- 
ness. From the dimuon and dielectron data, we obtain the upper limits on the cross section 
ratio Ra at 95% confidence level (C.L.) shown in Figs. Qupper) and (middle), respectively. 

In Fig. jlj the predicted cross section ratios for Zggj^ and Z'^ production are superimposed 
together with those for Gkk production with dimensionless graviton coupling to SM fields 
k/Mf\ = 0.05 and 0.1. The leading order cross section predictions for Zggj^ and Z'^ from 
PYTHIA using CTEQ6.1 PDFs are corrected for a mass dependent A:-factor obtained using 
ZWPRODP [34-37J to account for NNLO contributions. For the RS graviton model, a con- 
stant NLO /c-factor of 1.6 is used [38J. The uncertainties due to the QCD scale parameter and 
PDFs are indicated as a band. The NNLO prediction for the Z production cross section is 
0.97± 0.04 nb |126J. 

Propagating the above-mentioned uncertainties into the comparison of the experimental limits 
with the predicted cross section ratios, we exclude at 95% C.L. Z' masses as follows. From the 
dimuon only analysis, the Zggjyj can be excluded below 1027 GeV, the Z'^ below 792 GeV, and 
the RS Gkk below 778 (987) GeV for couplings of 0.05 (0.1). For the dielectron analysis, the 
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production of Zggjyj and Z'^ bosons is excluded for masses below 958 and 731 GeV, respectively. 
The corresponding lower limits on the mass for RS Gkk with couplings of 0.05 (0.10) are 729 
(931) GeV. 

7.1 Combined Limits on the Production Cross Section Using Dimuon and Di- 
electron Events 

The above statistical formalism is generalized to combine the results from the dimuon and 
dielectron channels, by defining the combined likelihood as the product of the likelihoods for 
the individual channels with forced to be the same value for both channels. The combined 
limit is shown in Fig. |4] (bottom). 

By combining the two charmels, the following 95% C.L. lower limits on the mass of a Z' res- 
onance are obtained: 1140 GeV for the Zggjyj, and 887 GeV for Z'^ models. RS Kaluza-Klein 
gravitons are excluded below 855 (1079) GeV for values of couplings 0.05 (0.10). Our observed 
limits are more restrictive than or comparable to those previously obtained via similar direct 
searches by the Tevatron experiments ||5l-0, or indirect searches by LEP-II experiments ||9l-[l2|, 
with the exception of Zggj^, where the value from LEP-II is the most restrictive. 

The distortion of the observed limits at ~400GeV visible in Fig.|4]is the result of a clustering 
of several dimuon and dielectron events around this mass. We have tested for the statistical 
significance of these excesses (p-values expressed as equivalent Z-values, i.e. effective number 
of Gaussian sigma in a one-sided test), using the techniques described in |39|. For the dimuon 
sample, the probability of an enhancement at least as large as that at 400 GeV occurring any- 
where above 200 GeV in the observed sample size corresponds to Z < 0.2; for the electron 
sample, it is less. For the combined data sample, the corresponding probability in a joint peak 
search is equivalent to Z = 1.1. 

In the narrow-width approximation, the cross section for the process pp — ?► Z' + X — > + X 
can be expressed IIT3l |34| in terms of the quantity c„zf ^ + c^Wrf, where c„ and contain the 
information from the model-dependent Z' couplings to fermions in the annihilation of charge 
2/3 and charge —1/3 quarks, respectively, and where Wu and contain the information about 
PDFs for the respective annihilation at a given Z' mass. 

The translation of the experimental limits into the {Cu,Cd) plane has been studied in the context 
of both the narrow-width and finite width approximations. The procedures have been shown 
to give the same results. In Fig. [s] the limits on the Z' mass are shown as lines in the {cd,Cu) 
plane intersected by curves from various models which specify {c^, Cu) as a function of a model 
mixing parameter. In this plane, the thin solid lines labeled by mass are iso-contours of cross 
section with constant Cu + (w^;/ Wu)c^, where W/i/ Wu is in the range 0.5-0.6 for the results rel- 
evant here. As this linear combination increases or decreases by an order of magnitude, the 
mass limits change by roughly 500 GeV. The point labeled SM corresponds to the Zggj^; it lies 
on the more general curve for the Generalized Sequential Standard Model (GSM) for which 
the generators of the LJ(1)t3j^ and U{1)q gauge groups are mixed with a mixing angle a. Then 
a = — 0.0727r corresponds to the Z'ggj^^ and a = and 7z/2 define the T^l and Q benchmarks, 
respectively, which have larger values of (c^, c„ ) and hence larger lower bounds on the masses. 
Also shown are contours for the Eg model (with x, ^, V, S, and N corresponding to angles 0, 
0.571, —0.29/1, 0.13/1, and 0.42/r, respectively) and Generalized LR models (with R, B — L, LR, 
and Y corresponding to angles 0, 0.5/r, — 0.13/r, and 0.25/r, respectively) [|34|| . 
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7 Limits on the Production Cross Section 




IVI [GeV] 

Figure 4: Upper limits as a function of resonance mass M, on the production ratio Ra of cross 
section times branching fraction into lepton pairs for Zggj^ and Gkk production and Z'^ boson 
production. The limits are shown from (top) the ji^ji^ final state, (middle) the ee final state 
and (bottom) the combined dilepton result. Shaded yellow and red bands correspond to the 
68% and 95% quantiles for the expected limits. The predicted cross section ratios are shown as 
bands, with widths indicating the theoretical uncertainties. 
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Figure 5: 95% C.L. lower limits on the Z' mass, represented by the thin continuous lines in the 
(cd, Cu) plane. Curves for three classes of model are shown. Colours on the curves correspond 
to different mixing angles of the generators defined in each model. For any point on a curve, 
the mass limit corresponding to that value of [c^, Cu) is given by the intersected contour. 

8 Summary 

The CMS Collaboration has searched for narrow resonances in the invariant mass spectrum of 
dimuon and dielectron final states in event samples corresponding to an integrated luminosity 
of 40 pb^^ and 35 pb^^, respectively. The spectra are consistent with standard model expec- 
tations and upper limits on the cross section times branching fraction for Z! into lepton pairs 
relative to standard model Z boson production have been set. Mass limits have been set on 
neutral gauge bosons Z' and RS Kaluza-Klein gravitons Gkk- A Z' with standard-model-like 
couplings can be excluded below 1140 GeV, the superstring-inspired Z'^ below 887 GeV, and RS 
Kaluza-Klein gravitons below 855 (1079) GeV for couplings of 0.05 (0.10), all at 95% C.L. The 
higher centre of mass energy used in this search, compared to that of previous experiments, 
has resulted in limits that are comparable to or exceed those previously published, despite the 
much lower integrated luminosity accumulated at the LHC thus far. 
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